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bX). We consider a simplicial minisuperspace model based on a cone over the triangulation 
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of the 3— sphere, in the presence of a massive scalar field, 0, with arbitrary scalar coupling 
term riRcf)^. By restricting all the interior edge lengths and all the boundary edge lengths 
to be equivalent and the scalar field to be homogenous on the 3— space, we obtain a two 
dimensional minisuperspace model {sj, 0j} for what is one of the most relevant triangulations 
of the spatial universe. We solve the two classical equations and find that there are both 
real Euclidean and Lorentzian classical solutions for any size of the boundary 3— space, 
After studying the analytic properties of the action in the space of complex edge lengths 
we then obtain steepest descents contours of constant imaginary action passing through 
Lorentzian classical geometries giving the dominant contribution and yielding a convergent 
wavefunction of the universe. We also show that the semiclassical wavefunctions for the 
Euclidean solutions associated with large boundary 3— spaces are exponentially suppressed. 
Consequently we can be confident that by using the SD contour associated with classical 
Lorentzian solutions the semiclassical approximation based on those classical solutions is 
well justified, clearly predicting classical spacetime in the late universe. This wavefunction 
is then evaluated numerically. 
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1 Introduction 



The sum over histories formulation of quantum gravity has undoubtably been one of the 
most useful tools in the long running program of quantising gravity. In this formulation 
an amplitude for a certain state of the universe is constructed by summing over a certain 
class of physically distinct histories that satisfy appropriate boundary conditions, weighted 
by their respective action. There are no a priori definitions of what the class of histories and 
boundary condition should be. As such what we have are proposals for both of them. In 
our opinion Hartle-Hawking's no boundary proposal formulated on [||] is the most natural 
one. As for the space of histories, until recently it was widely thought that these histories 
should be confined to smooth manifolds with well-behaved metrics. However, in and 0, 
Schleich and Witt put forward a powerful case for the generalisation of the space of histories 
to include smooth conifolds. Computations for concrete models in 0], and by us have 
reinforced that proposal. 

However, even this generalisation of the Hartle-Hawking proposal is plagued by at least 
two main problems. Firstly, the Euclidean gravitational action is not bounded from below 
which leads to the divergence of the Euclidean path integral. Secondly, there is no clear 
prescription for the correct integration contour to use. In Hartle proposed the use of the 
steepest descents contour in the space of complex metrics as the solution to both problems. 
Furthermore by choosing the steepest descents contour passing through the classical solutions 
of the theory, he made it very likely that the path integral be dominated by classical four- 
geometries, i.e., solutions of Einstein's equations and stationary points of the path integral, 
as desired for any wave function that is intended to represent our current Universe. Note that 
in this view, the fact that an integration solely over real-valued Euclidean geometries does 
not yield a convergent result for the path integral, is actually a good thing, for such a path 
integral would never predict the oscillatory behaviour in the late Universe that traditionally 
represents classical Lorentzian space-time. 

However, the usefulness of the above formulation in computational terms leaves some- 
thing to be desired. The resulting functional integrations over the metric tensor and matter 
fields are usually very hard, no matter what integration contour we choose. It is here that 
a simplicial formulation can be of great help. As we shall see below, in less than seven 
dimensions smooth manifolds are in a one-to-one correspondence with a special kind of sim- 
plicial complexes called combinatorial manifolds. This means that we can substitute a sum 



over smooth manifolds and metric tensors by a sum over simplicial complexes and their 
squared edge lengths, which in certain cases greatly simplifies calculations. However, as in 
the continuum case, a complete sum is still very hard, so we also end up using approximate 
minisuperspace models with features that greatly reduce the complexity of the calculations, 
in particular simplicial minisuperspace models based in Regge calculus. Such simplicial min- 
isuperspace models were introduced by Hartle . One typically takes the simplicial complex 
which models the topology of interest to be fixed and the square edge length assignments 
play the role of the metric degrees of freedom. The summation over edge lengths models the 
continuum integration over the metric tensor. This approach has several advantages. First 
by treating the four-geometry directly it is more adequate to deal with the Hartle-Hawking 
proposal, (with its four-dimensional nature), than the usual 3 + 1 ADM decomposition 
of space-time, where a careful study of how the four-geometry closes off at the beginning of 
the universe is essential. Second, by discretizing space-time the classical equations become 
algebraic which makes it easier to find classical solutions which are essential to the semi- 
classical approximation. Third, the simplicial minisuperpspace models offer the possibility 
of systematic improvement. 

In the present case we consider one such model where the topology has been fixed to that 
of a cone over 04, which is the simplest triangulation of the 3— sphere. As a matter sector 
we consider a massive scalar field with arbitrary scalar coupling to gravity, rjRcfP'. 

Unlike the results of our previous papers and |p where the only real classical solutions 
for universes with large ?> — D boundary, (like our own), were Lorentzian solutions, the intro- 
duction of the scalar coupling leads to the additional existence of real Euclidean solutions for 
universes with boundary of arbitrarily large size. However, the semiclassical wavefunction 
associated with such solutions proves to be exponentially suppressed and so the Lorentzian 
solutions still dominate, as they should if our model is to predict classical spacetime for the 
late universe. 

2 Simplicial Framework 

The crucial point in implementing any sum over histories formulation of quantum gravity 
is the specification of the set of histories to be considered. A history in quantum gravity is 
specified by its topology, smooth structure and geometry. In the case of the Hartle-Hawking 
approach, the histories considered have been topological spaces with the topology and smooth 



structure of a smooth compact manifold and a geometry specified on that manifold. Lately 
several papers like ||^ and have pointed to the advantages of extending such space of 
histories to include histories having more general topology, namely conifolds. In this paper 
we shall not deal with this issue. 

The traditional choice of histories described above derives from the fact that, in the 
(Euclideanized) classical theory of gravity a classical history is a Riemannian manifold 
(M", A,g), where g is a Riemannian metric and: 

Definition 2.1 The pair {M",A), where A = {{Ua,^a)} , is a smooth manifold with atlas 
A if it satisfies the following conditions: 

• Every point of has a neighbourhood Ua which is homeomorphic to an open subset 
of R^, via a mapping 

<^a:Ua^ R^ 

• Given any two neighbourhoods with nonempty intersection, then the mapping 

is a smooth mapping between subsets of R^. 

A topological space that satisfies only the first condition is called a topological manifold. 
Such spaces are not appropriate as histories because the lack of a smooth structure, i.e., 
atlas, makes it impossible to define essential concepts on them, like distance, continuous and 
differentiable functions (like scalar fields), integration, etc. 

The concrete implementation of a sum over smooth manifolds is very difficult. One of 
the main problems is how to provide a finite representation for the manifold-based histories. 
A simplicial formulation of quantum gravity aims to provide one such representation. For 
that to happen one must prove that there is a one to one correspondence between the set of 
smooth manifolds and some set of simplicial complexes. 

A simplicial complex somehow plays the role of topological manifold in the simplicial 
framework, in the sense that it also lacks the necessary structure to define essential concepts 
like dimension, distance, volume, curvature etc. 



Definition 2.2 A simplicial complex {K, \ K \) is a topological space \ K \ and a collection 
of simplices K, such that 

• \ K \ is a closed subset of some finite dimensional Euclidean space. 

• If a is a face of a simplex in K, then a is also contained in K. 

• If a a and ah are simplices in K, then a a Hai, is a face of both a a and a^. 

• The topological space \ K \ is the union of all simplices in K. 

If we are to be able to define essential concepts such as continuity and differentiability of 
functions on simplicial complexes we will need to introduce some kind of structure similar 
to that of smooth manifolds. To do so a few more definitions are necessary. 

Remember that the PL-join of a point a with a set L, denoted oL, is the union of all 
line segments joining points of L with the point a. This is also called a PL-cone over L with 
apex a. Remember also that: 

Definition 2.3 A map f : P ^ Q between two polyhedra P and Q, is said to be piecewise 
linear, (PL), if each point p in P has a cone neighbourhood N — aL such that 

f{Xa + iix) = A/(a) + nf{x) 

where x is in L, a is the apex of the cone N and A, /x>0, A-|-/x = l. 

Definition 2.4 A PL manifold M" is a topological manifold endowed with a PL-atlas A — 
{Ua, ^a)aeA> ^^^^ mapping 



is a piece-wise linear (PL) mapping between subsets of R\. 

These PL- manifolds are very closely connected with a special kind of simplicial complexes, 
the combinatorial manifolds, which are defined by imposing additional restrictions to the 
very general definition of simplicial complex. These restrictions make it possible to define 
concepts like distance, volume , curvature , etc. 
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Definition 2.5 A combinatorial n— manifold M.'^ , is an n— dimensional simplicial complex 
such that 



• It is pure. 

• It is n on- branching. 

• Any two n—simplices can be connected by a sequence of n—simplices, each intersecting 
along some {n — 1) — simplex. 

• The link of every vertex is a combinatorial (n — 1)— sphere. 

Indeed it can be shown that PL manifolds are equivalent to combinatorial manifolds, 
0]. So every combinatorial manifold admits a PL-atlas (f/a, $a)agA' ^"7 which concepts like 
continuity and differentiability of scalar fields in combinatorial manifolds can be defined, 
thus playing a similar role to that of the smooth structure in smooth manifolds. 

Furthermore given these definitions, following it can be proven that in less than seven 
dimensions every PL-manifold has a unique smoothing so we can state a very important 
result, namely: 

In less than seven dimensions, every smooth manifold, M", is triangulated by a unique 
combinatorial manifold, 

Obviously each smooth manifold has several distinct triangulations, what this result says 
is that only one of them is a combinatorial triangulation, i.e., a triangulation based on a 
combinatorial manifold, and not just any simplicial complex. 

So we see that the topological part of the sum over histories can be recast in terms of 
simplicial representatives of the "continuum" spaces. However in order to define a concrete 
sum over simplicial histories we still need to associate a metric and an action to each simplicial 
complex to be considered. Note that up to now we have not specified any kind of metric 
information associated with simplicial complexes. Once we have fixed the topology of the 
underlying simplicial complex, the most convenient way to attach metric information to it, 
is to use Regge calculus. 

2.1 Regge calculus 

A convenient way of defining an n— simplex is to specify the coordinates of its (n + 1) vertices, 
a = [0, 1,2, ...,n]. By specifying the squared values of the lengths of the edges Sij, we 
fix the simplicial metric on the simplex: 



g^^^^^^ = + (2) 

where i,j — 1, 2, ..n. 

So if we triangulate a smooth manifold M endowed with a metric g^i, by a homeomorphic 
simplicial manifold M , the metric information is transferred to the simplicial metric of that 
simplicial complex 

9i.Ax) — > 9ij{{sk}) = ^5l±^2 -hi (3) 

In the continuum framework the sum over metrics is implemented through a functional 
integral over the metric components {gn,^{x)}. In the simplicial framework the metric degrees 
of freedom are the squared edge lengths, and so the functional integral is replaced by a simple 
multiple integral over the values of the edge lengths. But not all edge lengths have equal 
standing. Only the ones associated with the interior of the simplicial complex get to be 
integrated over: 

I Dg,,{x) I D{s,} = n / Msi) (4) 

The boundary edge lengths remain after the sum over metrics and become the arguments 
of the wavefunction of the universe. In the simplicial framework the fact that the geometry 
of the complexes is completely fixed by the specification of the squared values of all edge 
lengths, means that all geometrical quantities, such as volumes and curvatures, can be 
expressed completely in terms of those edge lengths. 

We shall also be considering a scalar field with arbitrary mass and scalar curvature 
coupling, taking values 0^ at each vertex k of the complex. Generally these values will be 
independent, but like the edge lengths, not all have equal standing. Only the ones associated 
with interior vertices, {0^} are to be integrated over: 



Dm = llj d<f>i (5) 

The values of the field at the boundary vertices {4>f,} are just boundary conditions, becoming 
the arguments of the wavefunction. 

The Euchdeanized Einstein action for a smooth 4— manifold M with boundary dM, 
and endowed with a 4— metric, g^i,, and a scalar field with arbitrary mass m and scalar 
curvature coupling constant rj is 
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where K is the extrinsic curvature. 

Its simphcial analogue will be the Regge action for a combinatorial 4— manifold, M., with 
squared edge lengths {s^}, and with a scalar field taking values {(pv) for each vertex v of 
0]: 



—2 2A 



+ ^ E V^4(jV0? + ^ E V^4(j>i?,0? 

where: 

• (jfc denotes a fc— simplex belonging to the set of all fc— simplices in M.. 

• 0{&2)i is the deficit angle associated with the interior 2— simplex a\ = [u^] 



(a^)=27r- E ^d(4,^4) (6) 
o-4est(4) 
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and 9d{(j\, a^) is the dihedral angle between the 3— simplices era = [ijkl] and cig 
[ijkm], of (74 = [ijklm] that intersect at (72- Its full expression is given in Q. 

^/'(crg) is the deficit angle associated with the boundary 2— simplex a\: 



^(a^) = 7r- E ^^<i(^2,^4) (7) 

(T4e5t{(7^) 

VA:(crfc) for /c = 2,3,4 is the /c— volume associated with the /c— simplex, cxfc, and once 
again their explicit expressions in terms of the squared edge lengths are given in [HI. 



• V4(o"i), is the 4— volume in the simphcial complex M, associated with the edge ai, i.e., 
the volume of the space occupied by all points of Ai that are closer to (Xi than to any 
other edge of At. The same holds for V4(j) where j represents all vertices of Ai. 

• It can be shown, |Ty], that 

Um = ^ E V^2(^2)X(^2) 

cr2[jkl] 

where the sum is over all triangles a2 that contain the vertex j, and x{^2) is the deficit 
angle associated with (T2. We use the new symbol x because the triangle (J2 can be an 
interior or boundary triangle. 

It is easy to see that both V^i^ai) and V^iJ), can be expressed exclusively in terms of the 
edge lengths {s^}. In fact all the previous terms can be written exclusively in terms of {sk} 
and {0fc}. 

So we see that any history in QG of the type (M'*, A, ^f^,^, 0), where represents a 
topological manifold endowed with a smooth structure A, metric g^i,, and in the presence 
of matter fields represented by 0, has an unique simplicial analogue, (A1^, {sk}, {<Pj})- This 
allows us to concretely implement the formal sum over histories in terms of this finite rep- 
resentation as: 

where 

• {sj} are the squared lengths of the interior edges 

• {sh\ are the squared lengths of the boundary edges 

• {(pi} are the values of the field at the interior vertices 

• {0fc} are the values of the field at the boundary vertices 

Although the functional integral over metrics has been written explicitly in terms of the 
edge lengths, this expression is still heuristic because we still need to specify the list of 
suitable simplicial complexes Ai'^ we intend to sum over, the measure, and the integration 
contour to be used. To circumvent these problems we shall compute the sum approximately 



by singling out a subfamily of simplicial histories described by only a few parameters and 
carrying out the sum over these histories alone. 

An example of this is to adopt a simplicial minisuperspace approximation. We now 
describe in some detail the minisuperspace model we shall consider. 

3 Simplicial Minisuperspace 

We shall reduce our attention to a significant subfamily of simplicial histories characterised 
by the following restrictions: 

We shall consider that the universe has only one boundary and it is well approximated 
as a simplicial cone over the closed combinatorial 3-manifold a^, which is the simplest tri- 
angulation of the 3— sphere, S^. 

M^ = a*ai (9) 

The combinatorial manifold has been described in detail elsewhere p. We can see a 
representation of it in figure 1. It has 5 vertices, each connected to all others. 

Note that since all vertices of A^'^, (even the interior one) have combinatorial links that 
are homeomorphic to a 3— sphere, Ai'^ is a combinatorial 4— manifold. 

• By using a cone-like structure, translated by the existence of only one interior vertex, 
the apex a which we shall henceforth denote as 0, the only boundary of Ai"^ is . 
So it is very easy to define boundary simplices and interior simplices. If a simplex 
contains the interior vertex it is an interior simplex if not it is a boundary simplex. 
Moreover, all interior p— simplices in are cones over {p — 1)— simplices of ^4 with 
apex 0. 

• If we label the five boundary vertices of Ai'^ simply as 1,2,3,4,5, then the cone- 
like structure of leads to all interior edges being of the same form [0,6], with 
b = 1,2,3,4,5. So it makes sense to introduce the restriction that all interior edges 
have equal lengths whose squared value is denoted Si = sob- A similar assumption 
is made with respect to the boundary edge lengths, i.e., we consider them all to 
be equal to a common value Sij = Sb, with i,j = 1,2,3,4,5. We thus obtain an 
isotropic and homogeneous triangulation of the 4— universe. This leads to an enor- 
mous simplification in the metric part of the integration for the wavefunction (H), 



since the multiple integral / D{si} is reduced to a single integral / dsi. It also greatly 
simplifies the expression of the simplicial action since there will only be one type of 
boundary and interior triangle. 

• The simplifications assumed with respect to the edge lengths make it natural to assume 
that the scalar field is spatially homogeneous and isotropic. So we assume that the 
scalar field takes the same value (pb for all boundary vertices of A^^. The value at 
the interior vertex, 0j, is independent. Again, this leads to an simplification in the 
matter fields part of the integration for the wavefunction (|^), since the multiple integral 
/ D{(f)i}is reduced to a single integral / dcpi. 

3.1 Minisuperspace Wavefunction 

We can now concretely implement a simplicial minisuperspace approximation to the wave- 
function of the universe of the type (|^), as 



The Regge action for this minisuperspace can now be calculated. For simplicity we 
introduce rescaled metric variables: 



where H'^ = PA/3, and P = IQnG is the Planck length. We shall work in units where 
c = h = l. 

Then the volume of the 4— simplices in is 




(10) 
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(13) 



The volume of the 10 internal 2— simplices, in i 



IS 




(14) 



The volume of the 10 boundary 2— simplices, in i 



IS 



The volumes of the internal and boundary 3— simplices of M.^ are, respectively 



V2P 



(15) 
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(17) 
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There is only one kind of interior 2— simplex and boundary 2— simplex The dihedral angle 
associated with each interior 2— simplex is 



0{al) — arccos 



2^-1 
6^-2' 



(18) 



for the boundary 2— simplices we have 



^{^2) = arccos 



(19) 



With respect to the matter terms, the kinetic term vanishes when the edges (72 are 
boundary edges. The only non-vanishing contribution comes from the internal edges (72 = 
[Oj]. 

Computing the relevant volumes associated with the internal edges and all the vertices 
we conclude that the Regge action for this simplicial minisuperspace is 



^[C, S, (f)i, < 



+ 



TT — 2 arccos 



2y/W^ 



27r — 3 arccos 
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Note that we will be expressing the values of the field (f) and its mass m, in Planck units 

Thus in order to approximate the formal sum over histories by a fully computable ex- 
pression 



^[«4;s6;06]= / DiD<p,e-'^'^^-^^'^'t'^\ (20) 
Jc 

we only need to specify the integration contour C, and the measure of integration D^Dcpi. 

As in the previous cases studied by us p, 0, in this simphfied model the result yielded 
by a certain contour C is not very sensitive to the choice of measure if we stick to the usual 
measures, i.e., polynomials of the squared edge lengths. In our case we take 

ds- 9 

D(D4>. = ^<iA = ^df*. (21) 

Since in the case of closed cosmologies there is as yet no known explicit prescription for 
the integration contour, one usually takes a pragmatic view, in which we look for contours 



that lead to the desired features of the wavefunction of the universe. Following [jrT|, these 
features are: 

• It should yield a convergent path integral 

• The resulting wavefunction should predict classical spacetime in the late universe, i.e, 
oscillating behaviour when the is well approximated by the semiclassical approxi- 
mation. 

• The resulting wavefunction should obey the diffeomorphism constraints, in particular 
the Wheeler-DeWitt equation. 

It is well known that any integration contour over real metrics would yield a wavefunction 
that does not satisfy any of these basic requirements. On the other hand, an integration 
contour over complex metrics can, if wisely chosen, lead to a wavefunction that does satisfy 
them. 

In the simplicial framework, complex metrics arise from complex-valued squared edge 
lengths, (0). The boundary squared edge length, S, has to be real and positive for obvious 
physical reasons. But the interior squared edge length, ^, can be allowed to take complex 
values. 

Given the above requirements, 0] proposes that we use a steepest descents integration 
contour on the space of complex valued interior edge-lengths passing through the classical 
solutions that should dominate the wavefunction in the late universe. Before we test his 
proposal it is essential that we do the analytic study of the action as a multivalued function 
of the complex variable C,- 



3.2 Analytic Study of the Action 

The action is trivially analytic with respect to the variables 0i, 06 and S. But its dependence 
on the complex- valued ^ is much more complicated. So we shall investigate the analytic 
properties of / as if it was a function of ^ only, / = /[^], the other variables acting as 
parameters. 

The function has singularities at ,^ = and ^ = 1/3, and square root branch points 
at ^ = 1/4, 1/3 and 3/8. These branch points correspond respectively to the vanishing of 
the volume of the internal 2— simplices, 3— simplices and 4— simplices. Using 

arccos z — —i log{z + V z'^ — 1) 

we see that ^ = 1/3 is also a logarithmic branch point, near which the action behaves 
hke : 

I ~ i2 (5 V3 - lV3v4>ll') (1^) log (3e - 1) (22) 
The multivaluedness of associated with these branch points forces us to implement 
branch cuts in order to obtain a continuous function. In general, for terms of the type 

^yz — Zq we consider a branch cut (— c>o, 2:0]. So the branch cuts associated with the terms 
— 3/8, — 1/4 and ~ 1/3, altogether lead to a branch cut (— cxo, 3/8]. On the other 
hand, terms of the type arccos(^) have branch points at z = — 1,+1,(X), and usually the 
associated branch cuts are chosen as (—00, —1] U [l,+oo). These terms are also infinitely 
multivalued. 

The corresponding cuts for the term arccos are (|, |] U [|, |). On the other hand, 
associated with the term arccos — 7^= we have one cut §1 associated with arccositf^;), 

2^6^-2 V3' 8J V /' 

and another (—00, |] associated with u{z) — \^6^ — 2. 

So when we consider all these branch cuts simultaneously, we see that one way to obtain 
a continuous action / as a function of ^, is to consider a total branch cut (—00, |]. Note 
that this also takes care of the singularity at ^ = 0. Although the action then becomes 
a continuous function of ^ in the complex plane with a cut (— cxd, |], it is still infinitely 
multivalued. As usual in similar cases, in order to remove this multivaluedness we redefine 
the domain where the action is defined, from the complex plane to the Riemann surface 
associated with I. The infinite multivaluedness of the action reflects itself in I having an 
inflnite number of branches with different values. The Riemann surface is composed of an 
inflnite number of identical sheets, C — (— 00, |], one sheet for each branch of /. 



We define the first sheet Ci of I[^] as the sheet where the terms in arccos(2) assume 
their principal values. So the action in the first sheet will be formally equal to the original 
expression. Note that with the first sheet defined in this way, for real > 3/8 the volumes 
and deficit angles are all real, leading to a real Euclidean action for ^ G [|, +00) on the first 
sheet. On the other hand, when ^ is real and less than 1/4 in the first sheet , the volumes 
become pure imaginary and the Euclidean action becomes pure imaginary. For all other 
points of this first sheet the action is fully complex. 

When the action is continued in ^ once around all finite branch points (.^ = 1/4, 1/3, 3/8), 
we reach what shall be called the second sheet . It is easy to conclude that the action in this 
second sheet is just the negative of the action in the first sheet. 

Since by @ we see that the simplicial metric in each 4— simplex is real iff ^ is real, 
then the simplicial geometries built out of these 4— simplices will be real when ^ is real. 
Furthermore the corresponding eigenvalues of gij are A = {4(^ — 3/8), 1/2, 1/2, 1/2}, ^. 
So we see that for real C, > 3/8 we have real Euclidean signature geometries, with real 
Euclidean action, and for real ^ < 1/4, we have real Lorentzian signature geometries with 
pure imaginary Euclidean action. 

3.3 Asymptotic Behaviour of the Action 

If we are to compute the integration of along an SD contour, one of the essential things we 
have to know is the behaviour of the integrand, i.e., of the action, at infinity with respect to 
the variable ^. Only then can we be confident that the integral converges, with the classical 
solutions dominating the wavefunction for the late universe. 
It is easy to see that as ^ ^ 00 the action behaves like 
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(23) 



where 



S. 
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[27r — 3 arccos (1/3)] 



+ 
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(24) 
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The asymptotic behaviour of / for large ^ depends on whether or not S is larger than 
the critical value Sent- However, contrary to the corresponding model in [^] where there was 
no scalar curvature coupling and the critical value of S was 



_ 10[27r - 3arccos(l/3)] 

'-'crit — / \ 2 / \ y^^"^) 

. 1 



and as such could only take values in a limited range, now Scrit can be arbitrarily negative 
or positive because of its t] dependence. This means that the convergence of the integral 
along the SD contour, for a given S, will depend on the value of the coupling constant rj. 

4 Classical Solutions 

The classical simplicial geometries are the extrema of the Regge action we have obtained 
above. In our minisuperspace model there are two degrees of freedom ^, 0i, so the Regge 
equations of motion will be: 

dl , , 

= » (26) 



and 

dl 



0. (27) 



They are to be solved for the values of ^, 0,, subject to the fixed boundary data S*, 0b. 
The classical solutions will thus be of the form ^(S", 0;,), and 04(5", 0;,). The solution ^{S, (f)h) 
completely determines the simplicial geometry. 

Note that we shall be working on the first sheet. Of course, since on the second sheet 
the action is just the negative of this, the equations of motion are the same. And obviously 
every classical solution ^j{S, (pb) located on the first sheet will have a counterpart ^jj of the 
same numerical value, but located on the second sheet, and so with an action of opposite 
sign, I[^j(S,(j)b)] = (pb)]- So the classical solutions occur in pairs. 

The classical equation (pTl) is 
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A(0 = 1 + 60r/e A 
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(29) 



Introducing this equation into the first one (p^, we obtain a very long cubic equation in 
S for each value of ^, given fixed rj and 0^. 



where 



^3(0^' + A2iOS' + A^iOS + AoiO = 0, 



(30) 
(31) 



A2iO = -240\/2fl- 
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with = l/2{m/Hf. 

This equation can then be solved numerically for S", and by inverting the resulting so- 
lutions we obtain several branches of solutions ^ = ^ci{S,(f)b)- For obvious physical reasons 
we shall accept only solutions with real positive S. In figure 2 we show such solutions for 
7] = 0.015, m = 1, and (pt, = 1. Between the critical points ^ = 1/4 and ^ = 3/8 there are no 
real solutions. We chose a small value of rj and included negative branches so that the be- 
haviour of the solutions near the critical points is clear. The behaviour of the corresponding 
solutions for larger values of rj is of the same type, but the separation between the several 
branches is not so obvious. 

It is easy to see that for any positive value of S there is at least one classical Lorentzian 
solution (^ < 1/4). Furthermore, when S is large (late Universe) there is only one Lorentzian 



solution, C,i {S, (pb), (in the first sheet, of course), and is located near the critical point ,^ = 1/4, 
just like the results obtained in 0. There is also, of course , its counterpart in the second 
sheet ^^^(5", 0b), which though numerically equal has symmetrical action. However, the 
asymptotic behaviour of the classical solutions is very different from that in [0. Unlike in 
0, the Euclidean branch > 3/8), is not limited to a finite range of S. In the present case 
the Euclidean branch goes all the way to +oo, although this is not at all evident from figure 
2 because it is a large-scale behaviour, and this is a small-scale picture. We can see this 
large-scale behaviour in figure 3, but at the cost of the Lorentzian peak at ^ = 1/4 becoming 
indistinct from the imaginary ^ axis. In figure 3 we see that for each positive value of S 
there is always one pair of Euclidean signature solutions {S, (pb) = ^ii{S, (pb) £ [3/8, +oo), 
and ^ +00 as S" ^ +oo 

If we look in the opposite direction, i.e. C, — ^ — oo another surprise awaits us. The positive 
Lorentzian branch eventually becomes negative and connects with the upper negative branch. 

So the main difference introduced by the scalar coupling is the existence of Euclidean 
C,j{S,(f)b) = ^ii{S,(f)b) G [3/8, +oo) solutions for all positive values of the boundary edge 
length 5*. In a semiclassical analysis this seems to mean that the Lorentzian universe can 
nucleate with arbitrary size from an Euclidean regime. However the Regge action grows 
very fast as ^ goes from 3/8 to +oo, as we can see in figure 4. Consequently, the Euclidean 
classical solutions for large S are very strongly suppressed. 

For negative values of rj the solutions obtained have a very different behaviour, but they 
agree with the positive rj solutions in some very important points. In order to present these 
solutions more clearly we have separated the Lorentzian range < 3/8), presented in figure 
5, from the Euclidean range > 3/8), in figure 6. As can be seen in figure 5, where 
1] = —0.15, despite the somewhat bizarre behaviour, we still have that for large values of 5* 
there is an unique pair of classical Lorentzian solutions ^i{S, (pb) = ^ni^^ 'Pb) and these will 
be located near the critical point ^ = 1/4. Also, in the Euchdean regime, ]3/8,+oo), there 
is an unique pair of Euclidean solutions (S", cpb) = ^ii{S, (pb) £ [3/8, +oo) for any positive 
value of S, and those solutions go to +oo as S ^ -|-oo, as shown in figure 6. 

5 Steepest Descents Contour 

After studying the analytical and asymptotic properties of the action we can now focus on 
the Euclidean path integral that yields the wave function of the Universe (20), (21). 



As we have mentioned above there is as yet no universally accepted prescription for the 
integration contour C to use in quantum cosmology. Following Hartle we shall accept 
that the main criteria any contour should satisfy are that it should lead to a convergent path 
integral and to a wave function predicting classical Lorentzian spacetime in the late Universe. 
The steepest descents contour over complex metrics seems to be the leading candidate. In 
the simplicial framework, complex metrics arise from complex-valued squared edge lengths, 
(0). 

We shall look for the steepest descent (SD) contour, thinking of the action as a function 
of only the complex variable ^, and for the moment consider 0j to be only a real parameter in 
/ = /[^] to be integrated over later. In general, an SD contour associated with an extremum 
ends up either at oo, at a singular point of the integrand, or at another extremum with the 
same value of Jm(J). We have seen that when S is big enough the only classical solutions are 
a pair of real Lorentzian solutions (^/(>S', = ^^j{S,(f)h) — * 1/4), and a pair of Euclidean 
solutions (^^(S*, 0b) =?//(5', 0fc) >> 3/8). 

In both cases for each pair of solutions one of these solutions is located on the first sheet 
and the other on the second sheet. Thus, they have pure imaginary actions of opposite sign 
in the Lorentzian case and symmetric real actions in the Euclidean case. Given that their 
actions are different valued, in each pair there is no single SD path can go directly from one 
solution the other extremum, but it is possible for the two sections to meet in infinity, and 
together they form the total SD contour. In effect given that 

m = mr]]* 

and 

n^] = -nu 

where * denotes complex conjugation, we see that the SD path that passes through will 
be the complex conjugate of the SD path that passes through So the total SD contour 
will always be composed of two complex conjugate sections, each passing through one ex- 
tremum, and this together with the real analyticity of the action guarantees that the resulting 
wavefunction is real. 

For large S we choose to consider the SD contour associated with the Lorentzian solutions 
for two reasons. First, it is the only one likely to describe a late universe like our own. 
Second, the Euclidean action of the Euclidean solutions becomes very large very fast when 5" 



increases, which strongly suppresses these solutions in any computation of the wavefunction, 
except when S is small, (see figure 4). 

The SD path in the complex ^ Riemann surface of /, passing through a generic classical 
solution {C,ci{S, (j)b), 4>f{S, (pb)} is defined as: 

CsD{S,<PbA^) = {(e e R) : /m[J(^,e,0.,0.)] = i[US,<Pb),^f{S,^t)]] (32) 

where R is the Riemann sheet of the action, and /(^) = 

In figure 7 we show the result of a numerical computation of this path for m = 1, 06 = 1, 
T] = 0.015, S = 150, and (pi = 0.15. 

The SD path associated with the the other solution in the second sheet is just the mirror 
image of this, relative to the real ^ axis. Together they form the SD contour we are looking 
for. 

Changes in the values of 0b, 0j and rj, do not alter the generic shape of this SD contour.lt 
starts at +oo in the first quadrant and ends in the real ^ axis precisely at the classical 
solution to which it is associated. 

Going upward from the Lorentzian extremum on the first sheet, the SD contour proceeds 
to infinity in the first quadrant approximately along the parabola 
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The convergence of the integral along this part of the contour is dependent on the the 
asymptotic behaviour of the real part of the Euclidean action on the first quadrant of the 
first sheet 
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When there was no scalar curvature coupling, as in P], there was a finite maximum 
value that Scrit could take, no matter what the values of the parameters m and (ph were. 
This guaranteed convergence of the integral in the first quadrant, from a certain value of S", 



whatever the values of m and (pb- Now the situation is different because Scrit, (|2^), includes 
a term in rj that makes Scrit +00 as rj —00. 



6 Semiclassical Approximation 

One of the main requirements on any model is that it yields a wavefunction that in the 
late universe predicts a classical (Lorentzian) spacetime, like the one we experience. Now, a 
wavefunction of the universe will predict a classical spacetime where it is well approximated 
by the semiclassical approximation associated with Lorentzian classical solutions. From 
what we have seen above, the SD contour passing through the classical Lorentzian solutions 
satisfies that condition for large enough S. 

In our model there are two integration variables ^ and (pi and the full wavefunction of 
the universe is given by (JT^). We work under the assumption that (pi is to be integrated over 
real values, and ^ over the complex Riemann surface of the Euclidean action /. We now 
know that the integral over ^ can be calculated as a steepest descent (SD) integral for all 
the relevant values of (pi, and that the action peaks about the classical solutions ^. 

We can thus replace J^^^ d^e~^ by its semiclassical approximation based on the relevant 
classical extrema. This will give rise to Laplace type integrals in 0j, when the extrema 
have real Euclidean action, and to Fourier type integrals in 0j, when the extrema have pure 
imaginary Euclidean action. These integrals can then be shown to be dominated by the 
stationary points of the integrand which coincide with the classical solutions 0j, where 



This justifies the validity of the semiclassical approximation to the full wavefunction. 
So given the full wavefunction 



with an SD contour associated with real classical Lorentzian solutions {^^(5*, (p^)} with pure 
imaginary actions Ik = iI[^k{S,(ph);(pi\ = iIk{S, (pb, (pi), the semiclassical approximation of 
the wavefunction will be 



^m<p^] 

d(pi 








(33) 



where ' means derivative with respect to ^, and fik = sgn{l"). 



When the dominant extrema are real Euchdean solutions {^^{8, with real Euclidean 
actions, then after the semiclassical evaluation of the integral over ^, we are left with Laplace- 
type integrals over 0j which are dominated by the contributions coming from the stationary 
points of /[.^^(S*, 0b), which are precisely the classical solutions 0^(5', 0f,). So the semi- 
classical wavefunction will then be 



Since we are mainly interested in knowing if this model predicts classical Lorentzian 
spacetime for the late universe, we have computed the semiclassical wavefunction associated 
with the classical Lorentzian branch of solutions near ^ = 1/4 in figure 2. We considered 
7] = 0.225, m = 1 and 0f, = 1, and the result obtained is shown in figure 8. The behaviour 
exhibited during the late universe, (i.e. large values of S) is typical of that of a wavefunction 
describing a classical Lorentzian universe , as desired. 

As for the Euclidean solutions, their semiclassical contribution is exponentially suppressed 
except for small values of S; see figure 9. Furthermore, it should be noted that this suppres- 
sion becomes increasingly strong as r] grows. The peak in the semiclassical wavefunction is 
not caused by the behaviour of the Euclidean action, which is monotonically increasing as 
S increases, but by the pre-factor involving the second derivative of the action. 

It is clear that although there are classical Euclidean solutions for any value of 5* the 
probability associated with a Euclidean universe with large boundary, (large S), is very 
small, and so the late universe wavefunction should be well approximated by the semiclassical 
wavefunction associated with the Lorentzian solutions. 

7 Conclusions 

The addition of the arbitrary scalar coupling ?7i?0^ has produced results never seen in all of 
the previous simplicial models considered in [0, 0], and First, we note the existence 
of real Euclidean classical solutions for any size of the boundary t/iree— space. However, 
the contribution of these solutions to the wavefunction of the universe was seen to be ex- 
ponentially suppressed except for small boundary t/iree— spaces. In the late universe the 




(34) 



wavefunction was found to be dominated by the contribution of classical Lorentzian solu- 
tions. Second, the SD contours passing though these solutions are different from the ones 
obtained in 0. The SD path associated to the Lorentzian solution in the first sheet of 
Riemann surface of /[^], starts off at infinity moving downwards through the first quadrant, 
ending up precisely at the classical Lorentzian solution, and does not cross into the second 
sheet as before. Third, the behaviour of the real part of the action in the first quadrant now 
depends on the value of rj relative to the value of S", because Scrit +00 as 77 — — 00. So 
we see that since the value of r] is arbitrary so is the value Scrit from which the SD path 
converges. 

Nevertheless for any given 77 there is a value of S from which the SD integral is convergent 
and is dominated by the contribution of the Lorentzian classical solutions. Furthermore, 
larger values of S lead to a stronger peak in the action around those classical solutions 
and this makes the semiclassical approximations of the SD wavefunction quite good. The 
oscillatory behaviour of the semiclassical wavefunctions indicates that the wavefunction of 
the universe for this model predicts classical Lorentzian spacetime for the late universe, 
(large S). 
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Figure 2: Classical solutions ^{S, (/)b) for a scalar field of mass m = 1, (/){, = 1 
and scalar curvature coupling rj = 0.015. 
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Figure 3: Long-range view of the classical solution ^(5, 0^) in figure 2. 
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Figure 4: Euclidean action of the Euclidean classical solution ^(5, (f>b) for 
(f)b = l,m = 1, and scalar coupling rj = 0.225. 
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Figure 5: Lorentzian classical solutions ^{S > 0,(/)b), in the presence a scalar 
field of mass m = 1, (ph = 1 and negative scalar curvature coupling rj = —0.15. 
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Figure 6: Euclidean classical solutions ^{S > 0,^h), in the presence a scalar 
field of mass m = 1, (ph = 1 and negative scalar curvature coupling rj = —0.15. 
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Figure 7: Steepest descents contour when S = 150, (pb = 1,7] = 0.015, m = 1 
and = 0.15. Starting at infinity in the first sheet it ends at the real 
Lorentzian extremum at ^{S = 150, (f)b = 1) = 0.2222. Note that this kind of 
behaviour is maintained for other values of (pi. 
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Figure 8: Semiclassical wavefunction associated with a SD contour similar 
to the one in figure 7, where m = 1, (/){, = 1, but with r] = 0.225. We assume 
H = 7. 
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Figure 9: Semiclassical wavefunction associated with the Euchdean classical 
solutions, for a scalar field with m = 1, (/){, = 1, and scalar coupUng r] = 0.015, 
for H = 7. 
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